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We present exactly solvable examples that topological Mott insulators can emerge from topo- 
logically trivial states due to strong interactions between atoms for atomic mixtures trapped in 
one-dimensional optical superlattice systems. The topological Mott insulating state is characterized 
by nonzero Chern number and appears in the strongly interacting limit as long as the total band 
filling factor is an integer, which is not sensitive to the filling of each component. The topological 
nature of the Mott phase can be revealed by observing the density profile of the trapped system. 
Our results can be also generalized to the multi-component atomic systems. 
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Introduction.- Topological band insulators have at- 
tracted great attention in condensed matter physics since 
they were discovered in spin-orbit coupled two-and three- 
dimensional systems P, Ql ■ As the topological insulator 
(TI) is expected to be robust against weak perturbations, 
it is important to consider the effect of interactions on the 
topological band insulators Although generally a 

strong interaction may open the energy gap and break 
the TI, there exists a class of topological insulator called 
as the topological Mott insulator, where the interaction 
effects are responsible for TI behavior Q. The topo- 
logical Mott phases have attracted increasing attention 
|8l4ld| since the concept is proposed as interaction plays a 
crucial role in the formation of both the Mott phase and 
nontrivial topology. As almost all these results based on 
the mean-field approximation, examples with the inter- 
action effect counted exactly is particularly important for 
our understanding of the topological Mott insulator. 

In this work, we explore the realization of topologi- 
cal Mott phases in cold atomic systems trapped in one- 
dimensional (ID) quasi-periodic optical lattices, which 
can be generated by superimposing two ID optical lat- 
tices with commensurate or incommensurate wavelengths 
13|. Cold atomic systems in ID quasi-periodic lat- 
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tices have been extensively studied 
on the Anderson localization 12 1. 



17[ with a focus 
However, their non- 
trivial topological features are recognized only very re- 
cently [18j, JJi] . Particularly, with the experimental ob- 
servation of the topological edge states in ID photonic 
quasi-crystals fl9'| there is a growing interest in the study 
of topological properties in the ID quasi-periodic lattices 
[20|-|25|. It has been shown that the free fermion system 
with its sub-bands being fully filled is a topological non- 
trivial insulator, which can be characterized by a nonzero 
Chern number in the extended two-dimensional (2D) pa- 
rameter space [Til . In this work, we study the interacting 
multi-component atomic mixtures in the ID superlattice 
with its sub-bands are partially filled by atoms. In the 
absence inter-component interaction, the system may be 



a Fermi metal, a Bose superfluid or their mixture, de- 
pending on the component of mixture being fermion or 
boson. We find that a topological Mott insulator may 
emerge in the strongly interacting regime if the total 
band filling is an integer. This conclusion is exact in 
the infinitely strong interaction limit as it is based on an 
exact mapping which relates the many-body wavefunc- 
tion of two-component mixture to the wavefunction of 
free fermion system. Our study provides a simple way 
to realize the topological Mott insulator in cold atomic 
systems and may deepen our understanding of the topo- 
logical Mott insulator. 

Interacting atomic mixtures in optical superlattices. - 
We consider the ID atomic mixtures loaded in a bichro- 
matic optical lattice ll-l3|, which is described by i7 = 
Ha + Hi with 

Ho = -t {c\„h+\.a + H.c.) + ^ K;n,,<, (1) 



and 



where Vi — Acos(27raj-|-(5) with A controlling the strength 
of commensurate potential, a. tuning the modulation pe- 
riod and b being an arbitrary phase, Ci^a are bosonic or 
fermionic annihilation operators localized on site i, and 
"■io- = cj Ci^ . Here a 4- denotes the pseudo-spin in- 
dex of two components of atomic mixtures, which can 
be either fermion or boson. The interaction parame- 
ter U denotes the inter-component interaction strength 
and Ua denotes the intra-component interaction strength 
between bosonic atoms. For the Fermi-Fermi mixture, 
Ua — Q due to the Pauli exclusion. For the Fermi-Bose 
mixture, one of is zero if a denotes the fermion com- 
ponent. In principle, both V and Ur, can be tuned ex- 
perimentally by the Feshbach resonance. The hopping 
amplitude t is set to be the unit of the energy (t = 1). 
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In the absence of interaction, i.e., Hi = 0, the eigen- 
value equation, namely the Harper equation, is given by 

-[0„(i-|-l) + 0„(i-l)] + Acos(27rm + (5)0„(i) = e„(/)„(i), 

(3) 

where denotes the single particle wave function and 
e„ the n-th single particle eigenenergy. For a superlattice 
with a — p/q {p and q are co-prime integers), the system 
has a unit cell of q sites and the single-particle spectrum 
is split into q bands, which also changes periodically as 
the phase S varies from to 27r. Given that the number 
of CT-component atoms is Na and the number of lattice 
sites is L, we define the component-depending band fill- 
ing factor i^g. as = N^/Nceii with Nceii = L/q being 
the number of primitive cells. For a Fermi system, the 
system with the band filling factor v„ — m (m is an inte- 
ger smaller than q) corresponds to a band insulator with 
the lowest m bands being fully filled by the cr-component 
fermion. Such a band insulator has been demonstrated to 
be characterized by a nontrivial topological Chern num- 
ber in a 2D parameter space spanned by momentum and 
the phase of 5 jl^. 

In the following calculation we shall consider the two- 
component system with fractional filling factors and 
but the total band filling factor v = + being 
an integer. For a free two-component Fermi system with 
such a filling, the sub-band is only partially filled and the 
system is a topologically trivial conductor. We shall show 
that a Mott phase can emerge from the conducting phase 
with an integer total filling factor when the interaction 
effect is considered. 

Emergence of Mott phase for the Fermi-Fermi 
mixture.- To give a concrete example, we first consider 
the Fermi Hubbard model in the superlattice, i.e., the 
Hamiltonian ([T|) and ([2]) with Ua = 0. We consider the 
case with a — 1/i and v = 1, for which n = N/L = 
1/3. We calculate the charge gap, which is defined as 
A = [EoiN + 1) + Eo{N - l)]/2 - Eo{N), by numerically 
diagonalizing the Hamiltonian, where Eo{N) represent 
the ground state for the N-atom system. In the Fig. la, 
we show the charge gap versus the repulsive interaction 
strength U for the equal- mixing case with ly-^ = i^i = 1/2. 
It is shown that the charge gap increases with increasing 
U and tends to a constant in the large U limit. For sys- 
tems with different sizes, we find the charge gap in the 
large U limit approaches to the half size of the band gap 
Ab between the lowest band and the second one. The nu- 
merical results suggest that a Mott insulator is emergent 
in the strongly interacting limit as the gap is induced 
by the repulsive interaction between fermions of different 
components. 

Next we calculate the charge gap versus the repul- 
sive interaction strength U for various U' for the equal- 
mixing Fermi-Bose mixture. Here we use U' to repre- 
sent intra-component boson-boson interaction strength 
for the two-component Fermi-Bose mixture. As shown 
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FIG. 1: (Color online) The charge gap versus the interac- 
tion strength U for equal-mixing Fermi-Fermi system (a) and 
Fermi-Bose system (b) with a = 1/3, A = 1.5, 5 — 0, 
N-f — N_i = 2 and L = 12 under periodic boundary condi- 
tions (PBC). 



in Fig. lb, with the increase in the boson-boson inter- 
action strength, the charge gap increases and the curve 
approaches the curve of Fermi-Fermi mixture. In the 
limit of hard-core boson, i.e., [7^ — > oo for the Bose-Bose 
mixture and J7' — ^ oo for the Fermi-Bose mixture, both 
the Bose-Bose Hubbard model and Bose-Fermi Hubbard 
model can be mapped into the Fermi Hubbard model 
by extended Jordan- Wigner transformations [26^ . There- 
fore, the above discussion on the Fermi-Fermi mixture 
can be directly applied to hard-core-boson limit of the 
Bosc-Bose and Fermi-Bose mixtures. For these systems 
with v = 1, the Mott phase is emergent with increasing 
the strength of inter-specie repulsion U. 

Atomic mixtures in the strongly interacting limit- The 
emergence of the charge gap for the superlattice system 
with integer band filling factor can be clearly understood 
in the strongly interacting limit. In the infinitely strong 
repulsion limit, U — >■ oo, we can construct the many-body 
wavefunction of the interacting fermion system exactly 
by using the hard-core contact boundary condition [27| 
and group theoretical methods [2^ . Combining with the 
Pauli exclusion principle, the effect of an infinitely strong 
interaction can be reduced to a hard-core contact bound- 
ary condition (a;i , ai ; ■ • • \ XN,cfN) \xi=xj= 0, which 
does not depend on spin configurations. According to 
Ref. 28|, the many-body wave function of the system 
can be represented as ^ = ^^^'s, where the spatial 
wave function a is composed of Slater determinant of 
N = + orbitals 0i(x), • • • ,(j)N[x), given by 



iI)a{xi,...,xn)^{N\) ^det[ 



(4) 
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FIG. 2: (Color online) Schematic diagram for the ground state 
energy level occupations for system with a — 1/3, — Ni — 
2 and L — 12. (a) is for U = and (b)is for U oo. 



with (f)n (x) (x — ia with a the lattice constant) the eigen- 
state of the single particle Hamiltonian Ho, whereas ^5 
is a mapping function composed of linear combination 
of production of sign fmiction and basis tensor of spin, 
whose explicit form is given in Ref. [28]. Effectively, the 
effect of infinite repulsion is to generate a Pauli exclu- 
sion between different components of fermions. Such an 
effect has been recently observed by the experiment in 
a two-particle system with tunable interaction using two 
fermionic ^Li atoms [2^. The ground state of the system 
in the limit of t/ ^ 00 is composed of states with the 
lowest N single particle states occupied, and the ground 
energy is given by Eq = J2n=i ^ comparison, 

the ground energy of the free fermion system is given by 

Now we can understand the formation of Mott phase 
in optical superlattices by considering the limit case with 
infinitely strong interaction. As illustrated in Fig. 2, for 
the superlattice with a ~ 1/3, the energy levels split 
into three bands. When U = 0, the lowest band is only 
partially filled for 1^^ = = 1/2. However the lowest 
band is fully filled in the limit of f/ — >• cx) corresponding 
to = 1 as each energy level can only be occupied by 
a fermion with either spin up or spin down, and conse- 
quently a charge gap is opened. According to the defini- 
tion of charge gap, the charge gap equals to half of the 
band gap as adding a fermion costs the energy of Afc. 
This is consistent with the numerical result displayed in 
Fig.l. For a general case with a — 1/q, the energy lev- 
els are composed of q bands. As long as = m with 
m — 1, ■ ■ ■ ,q — I, a Mott insulator is formed as U ^ 00. 
For a fractional v, for example, v < 1, the lowest band is 
not fully filled even in the strong interaction limit, and 
thus no finite charge gap is opened. 

The above discussion can be directly applied to the 
Bose-Bose mixture and Fermi-Bose mixture with in- 
finitely strong inter-component and intra-component in- 
teraction. As C/ — >■ 00 and Uc — >■ 00, the effect of the 
infinitely strong interaction can be also reduced to the 
hard-core contact boundary condition ^'(2;^ = Xj) = 0, 
which enforces an effective Pauli exclusion to hard-core 
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FIG. 3: (Color online) Low-energy spectrum versus 5 for two- 
component mixture systems with A = 1.5 and different a in 
the infinitely repulsive limit, (a) a = 1/5, A*' — 30, L = 90 
under PBC; (b) a = 1/3, iV = 30, L = 91 under OBC; (c) 
a = 1/5, iV = 40, L = 100 under PBC; (d) a = 1/5, TV = 40, 
L = 101 under OBC. 

bosons and between different components of atoms. Con- 
sequently, the exact wave function of the system can be 
also represented as ^' = ^'a^'s, where 'if a is identical 
to the expression of ^ , but ^ 9 has different form for 
different kind of mixtures [l^, O, [s^ . For the Bose-Bose 
mixture, = IlKi j<Af *5"^(^« ~ where sgn(x) is 
the sign function. The explicit form of the mapping func- 
tion s for the Fermi-Bose mixture is given in Ref. (27j . 

Topological feature of the Mott phase.- To character- 
ize the topological feature of the many-body states, it is 
convenient to introduce a generalized boundary phase 9 
by applying the twist boundary condition along the ID 
chain. In the 2D parameter space of {9,5), we can cal- 
culate the Chern number of the many-body state, which 
is defined as an integral invariant C = 57 / d9dSF{6, 6), 
where Fi9J) = Im((||||f ) - (||||f » is the Berry 
curvature [3ll |32|. For the two-component mixtures in 
the strongly interacting limit, we notice that the many- 
body wave functions under the twist boundary condition 
can be represented as ^{9,S) = ^'^(0, (5)\l/s, where only 
^I'a varies with the change of 9 and 5, whereas is 
independent of 9 and S as it is composed of combina- 
tion of production of sign functions. This greatly sim- 
plify the calculation of the Chern number as F{9, 6) — 
M(^l^) - (^1^)), where is given by 
Eq.(|4l). Consequently, in the strong interacting limit, the 
Chern number for atomic mixtures composed of Nf and 
N_i two-component atoms is identical to the Chern num- 
ber of the system composed of N free fermions. For the 
strongly interacting atomic mixtures on the superlattice 
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with a = 1/3, the Mott insulators are formed when v = 1 
and = 2, corresponding to states with Chern number 
C = 1 and C = —1. These Mott insulators are thus 
topologically nontrivial insulators. 

According to the bulk-edge correspondence in general 
topological insulators, one may expect that the topologi- 
cal Mott insulators should display nontrivial edge states 
for the system with open boundary conditions (OBC). 
To see it clearly, we display the excitation spectrum for 
mixture systems in the infinitely strong repulsion limit 
with 1/ = 1, a = 1/3 and v = 2, a — 1/5 in Figl3]versus 
the phase S. The system with i/ = I and a = 1/3 is 
a topological Mott insulator characterized by the Chern 
number C — 1, whereas the topological Mott insulator 
corresponding to v — 2 and a = 1/5 is characterized by 
the Chern number C = 2. As shown in Fig. and Fig. 
[3};, there is an obvious gap between the ground state and 
the first excited state for systems with PBC. However, 
as shown in Fig. [Sb and Fig. |3ji, edge states appear 
in the gap regimes for systems with OBC. As the phase 
varies from — tt to tt, the edge states connect the ground 
state to the excited band. When the interaction deviates 
from the infinite limit, despite of the exact mapping no 
longer holding true, our numerical results indicate that 
the topological Mott phase still exits. As shown in FigHJ 
the excitation spectrum for the system with U = 100 
exhibits almost the same behavior as shown in Fig. 3b. 
Even when the interaction strength decreases to U = 5, 
the spectrum still has similar structure. The main dif- 
ference is that the ground state levels are broadened due 
to spin fluctuations and the degeneracy of ground states 
in the infinite interaction limit is lifted for finite interac- 
tions. 

Our results can be directly extended to the general 
multi-component atomic systems trapped in ID superlat- 
tices. When the inter-specie and inner-specie interactions 
between atoms go to the strongly interacting limit, an ef- 
fective Pauli exclusion between atoms arises. Therefore, 
a topological Mott insulator is expected to appear in the 
topologically nontrivial superlattice as long as the total 
filling factor = is an integer, where k is the com- 

ponent index of the multi-component atomic mixture. 

Experimental detection.- A possible way to observe 
the topological Mott insulator is to detect the density 
distribution of the superlattice systems with an exter- 
nal confining potential, i.e., Vi in Eq.(Il]) is replaced by 
Vi = Xcos{2nai + S) + Vnii — io)^, where Vh is the 
strength of the harmonic trap with being the position 
of trap center. In Fig. 5, we show the local average density 
distribution of two-component fermions with imbalance 
populations trapped in the optical superlattice with a 
harmonic trap. The local average density is defines as 
Pi = E*l-Af "*+i/(2^ + 1) with M < L [181, where 
Tii = {fii,-\) + {fi'i,i)- In the strongly interacting limit, two 
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FIG. 4: (Color online) Low-energy spectrum versus 5 for the 
Fermi- Fermi mixture with a = 1/3, A = 1.5 A^^^ = — 2, 
L — 13 and various U under open boundary condition, (a) 
U = l, {h)U = 5, (c) U = 10 and (d) U = 100. 
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FIG. 5: (Color online) The local average density distributions 
for the superlattice with a = 1/3 (a), and a — 1/5 (b) trapped 
in the harmonic trap. The system is with 1600 sites, Nf = 
400, A^j. = 500, A = 1.5, 5 = and Vh = 2 x 10"''. The dot 
line is for the case of [/ = and the solid line is for the case 
of f/ oo. Here, we take M = (g - l)/2. 



for a — 1/5. We note that locations of plateaus appear 
at p{a) — a,l — a, 2a, 1 — 2a, 4a, 1 — 4a, if the values 
are in the range of (0, 1). For contrast, the density pro- 
files for [/ = do not exhibit these plateaus. The Chern 
number can thus be deduced from the plateau distribu- 
tion by using the Streda formula [H, lH C = It 
is straightforward to get C — 1,-1 for pi a, 1 — a and 
C = 2, -2 for p, = 2a, 1 - 2a. 
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Summary.- In summary, we studied the realization of 
topological Mott insulators in interacting atomic mix- 
tures trapped in ID optical superlattices. We give a clear 
interpretation for the formation of topological Mott insu- 
lator by using an exact mapping, which relates the many- 
body wavefunction of atomic mixtures in the strongly in- 
teracting limit to that of the free fermion system. The 
topological Mott insulator displays nontrivial edge states 
and can be characterized by a nonzero Chern number. 
The topological feature of the Mott insulator can be re- 
vealed by detecting plateaus of density profiles of the 
trapped lattice systems. Our results pave the way for ex- 
perimentally studying topological Mott insulators in ID 
optical lattices. 
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